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ABSTRACT 

From Hubble Space Telescope images with 0'.'05 resolution we identify four stars brighter than V = 25 
mag within 2'.'5 of SN 1993J in M81 which contaminated previous ground-based brightness estimates 
for the supernova progenitor. Correcting for the contamination, we find that the energy distribution 
of the progenitor is consistent with that of an early K-type supergiant star with My ~ —7.0 ± 0.4 
mag and an initial mass of 13-22 Mq. The brightnesses of the nearby stars are sufficient to account 
for the excess blue light seen from the ground in pre-explosion observations. Therefore, the SN 1993J 
progenitor did not necessarily have a blue companion, although by 2001, fainter blue stars are seen in 
close proximity to the supernova. These observations do not strongly limit the mass of a hypothetical 
companion. A blue dwarf star with a mass up to 30 Mq could have been orbiting the progenitor 
without being detected in the ground-based images. Explosion models and observations show that the 
SN 1993J progenitor had a helium-rich envelope. To test whether the helium abundance could influence 
the energy distribution of the progenitor, we calculated model supergiant atmospheres with a range of 
plausible helium abundances. The models show that the pre-supernova colors are not strongly affected 
by the helium abundance longward of 4000 A, and abundances ranging between solar and 90% helium 
(by number) are all consistent with the observations. 

Subject headings: supernovae: general — supernovae: individual (SN 1993J) — galaxies: individual (M81, 
NGC 3031)— stars: evolution 



1. INTRODUCTION 

Direct identification of the progenitor star has been 
scarce for the more than 2000 known historical super- 
novae (SNe); hence, our knowledge of the types of stars 
that become SNe is based primarily on models of stellar 
evolution and observations of the environments in which 
SNe occur. Pre-SN stars have been detected for only 
five events: SN 1987A in the LMC (e.g., Walborn et 
al. 1987), SN 1961V in NGC 1058 (Goodrich et al. 1989; 
Filippenko et al. 1995), SN 1978K in NGC 1313 (Ry- 
der et al. 1993), SN 1997bs in NGC 3627 (Van Dyk et 
al. 1999), and SN 1993 J in M81 (Fihppenko 1993a; Alder- 
ing, Humphreys, & Richmond 1994, hereafter AHR; Co- 
hen, Darling, & Porter 1995). Interestingly, the proper- 
ties of all five of these SNe are atypical of type II super- 
novae (SNe II). In fact, the highly unusual properties of SN 
1961V and SN 1997bs (Van Dyk et al. 2000) suggest that 



they were actually eruptions of extremely massive stars, 
similar to rj Carinae, rather than genuine SNe. 

SN 1993J was initially classified as a SN II, because its 
optical spectra showed hydrogen lines (Filippenko 1993b; 
Garnavich & Ann 1993), but it soon developed the char- 
acteristics of a SN lb (Filippenko & Matheson 1993; Fil- 
ippenko, Matheson, & Ho 1993; Filippenko, Matheson, & 
Barth 1994) — hence, the classification as Type lib. Sev- 
eral groups (see the review by Wheeler & Filippenko 1996) 
concluded that the progenitor of SN 1993J was a massive 
star (10-20 Mq) that had lost most of its hydrogen en- 
velope before exploding. The models generally employ a 
close companion star to strip the progenitor of its envelope, 
since mass loss through a wind is not efficient at such low 
initial masses. The companion gains mass and may reach 
more than 15 Mq before the supernova explosion (Woosley 
et al. 1994). Hofiich, Langer, & Duschinger (1993) offered 
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an alternative model, in which a single, very massive star 
loses its envelope through a wind before exploding; the 
results, however, do not fit the observations as well as do 
binary star models. 

The photometric properties of the SN 1993 J progenitor 
were estimated by AHR, using an extensive ground-based 
photographic and CCD archive for M81. Their analysis 
was affected by the range in quality of the observations 
and the crowded field near the SN. AHR concluded that 
the progenitor was a KO supergiant with contamination 
from one or more blue stars projected near the SN and 
possibly light from an OB association that spawned SN 
1993J. 

From a set of CCD images obtained in the 1980s to 
search for novae in M81, Cohen et al. (1995) deter- 
mined that the SN 1993J progenitor was not variable at 
the 0.2 mag level. They also confirmed the AHR esti- 
mate of the progenitor color. Crotts (1995) obtained high- 
quality CCD images 600 days after maximum brightness 
and found a red star with / ^ 22.7 mag within 1" of 
SN 1993J. This star was responsible for the elliptical ap- 
pearance of the progenitor in the AHR set of images, even 
under good seeing conditions. 

In the course of the Supernova Intensive Study (SINS) 
collaboration, multiwavelength Hubble Space Telescope 
{HST) images had been obtained of SN 1993 J and the 
surrounding star field in 1994 and 1995, while the SN was 
fading. The SN and the field have also been observed in 
2001, though not as the main target, by program GO-9073 
(see Liu, Bregman, & Seitzer 2002). In this paper, we use 
all these data to help resolve some questions that were 
raised by the ground-based data about the progenitor. In 
particular, we estimate the contamination of the progeni- 
tor light by stars which are unresolved in the ground-based 
images and attempt to determine the source of the excess 
blue light seen by AHR. We also estimate the brightness, 
color, age, and initial mass of the progenitor. In §2 we 
describe the observations and photometric reduction. The 
implications of the HST data for the pre-SN observations 
are discussed in §3, and conclusions are presented in §4. 

2. OBSERVATIONS AND DATA REDUCTION 

Short exposures were obtained, when the SN was still 
fairly bright, on 1994 April 18 (UT is used throughout 

this paper) with the HST WFPC2 through fihers F336W, 
F439W, F555W, F675W, and F814W, which approximate 
standard Johnson-Cousins UBVRI bandpasses. After the 
SN faded further, images were obtained on 1995 January 
31 through the same filters, and additionally, the UV fil- 
ter F255W. Exposures through fihers F450W, F555W, and 
F814W were also obtained by GO-9073 on 2001 June 4. 
Table 1 summarizes all of the HST observations. 

Photometry of these images was performed using HST- 
phot version 1.1 (Dolphin 2000a,b). This package au- 
tomatically accounts for WFPC2 point-spread fimction 
(PSF) variations and charge-transfer effects across the 
chips, zero-points, aperture corrections, etc., and returns 
UBVRI magnitudes as output. HSTphot appears to have 
a limiting (image-to-image) photometric accuracy of ~0.02 
mag (Dolphin 2000a). The transformations, even at faint 
magnitudes, should be at this level or better (Dolphin 
2000b). 



Approximately 25 stellar sources are within 2'.'5 of the 
SN that can be identified in the sum of all images ob- 
tained in 1995 (see Figure 1). Of these, only six have 
B < 25 mag, and only four have V < 25 mag. For iden- 
tification purposes, the angular distances and position an- 
gles from SN 1993J for stars A F have been listed in Ta- 
ble 2. The magnitudes from each epoch, along with the 
adopted magnitudes for the stars in the SN environment 
(the uncertainty- weighted average of the three epochs), 
are given in Table 2. Measurements for the brightness 
of SN 1993J itself are given in Table 3. 

A distance to M81 of 3.63 ±0.34 Mpc {fi = 27.80 ±0.20 
mag) has been determined from HST observations of 
Cepheid variables (Preedman et al. 1994). The redden- 
ing to SN 1993J is not well-determined. Richmond et 
al. (1994) estimate E{B - V) = 0.08-0.32 mag. From 
early-time spectra, Clocchiatti et al. (1995) find a visual 
extinction Ay = 0.74 ± 0.05 mag. As AHR point out, 
the range in likely Ay to the SN is about 0-1.5 mag. We 
adopt Ay = 0.75 mag, in agreement with Clocchiatti et 
al. (1995) and consistent with the midrange of colors and 
magnitudes of stars in the SN's environment (see below). 
Therefore, a star with B = 25 mag, if physically associated 
with M81, would have absolute magnitude Mb ~ —3.6. 

The SN was still bright in 1995 January {V w 18.8 mag) 
and faint stars closer than ~0'.'25 (~4.4 pc at the distance 
of M81) would not be resolved, if present. Stars with 
B > 25 mag, even in 2001, are more than 11 times fainter 
in B than the estimated brightness of the progenitor (see 
below) and could not substantially affect the progenitor 
photometry. 

The position of the red star reported by Crotts (1995; his 
estimates are 69 = 0'.'84, P.A. = 344°) is consistent with 
Star A in the HST images. However, Crotts' measurement 
for the star, I w 22.7 mag, is ~ 2 mag fainter than our 
HST measurement of 20.7 mag. Also, the color for this 
star, R — I = 1.14 mag, based on the HST measurements, 
is much redder than the R — I « 0.5 mag reported by 
Crotts. Star A may be somewhat variable, as seen from the 
magnitudes derived from all three epochs of HST observa- 
tions. On the other hand, the variability is certainly not 
of this degree, nor can we easily explain the difference due 
to bandpass differences. This discrepancy demonstrates 
the difficulty in obtaining accurate photometry of closely 
blended stars, even from the best ground-based sites. 



3. DISCUSSION 

Analysis of archival ground-based images by AHR 

showed that no single star could fit the progenitor's energy 
distribution. They suggested that the spectrum is best fit 
by a slightly reddened KO supergiant, plus light from an 
OB association surrounding the SN. From the HST im- 
ages, we find four stars, A-D, with V <^ 24: mag, that 
fall within the resolution element of observations typical 
for AHR (FWHM w 1'.'5, corresponding to FWZI « 2'.'5), 
which may have influenced flux measurements of the pro- 
genitor by AHR. Star A appears to be variable in all bands, 
which might possibly have influenced the progenitor flux 
measurements by Cohen et al. (1995). However, they rule 
out variability for the progenitor at the ~0.2 mag level; the 
variability we measure in V for Star A is comparable to, or 
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slightly larger than, this level, although the measured be- 
havior here is consistent with that shown by Cohen et al. 
(We note that the apparent variability for Star A we mea- 
sure in B and -R is > 0.2 mag.) Two blue stars are seen 
in Figure 1 very close to the SN, but they are near the 
limit of the detection threshold and are therefore likely 
too faint to substantially affect the progenitor photome- 
try. The V = 23.5 mag star just east of due north and the 
V = 24.3 mag, somewhat extended, source to the south- 
west, both near the edge of the circle in the figure, are well 
in the wings of the typical ground-based seeing disk. All 
other stars within the circle are quite faint {V » 25-26 
mag). 

Based on a visual inspection of the HST images, SN 
1993J does not appear to be a member of a cluster or 
association. Using IRAF^/DAOFIND, a detailed compar- 
ison between the density of stars within 2" of the SN with 
the overall star density on the B-band frames also does not 
indicate a statistically significant stellar clustering around 
the supernova. SN 1993J is at a large radial distance from 
the center of M81 (4.'8), located along an outer spiral arm, 
which is interspersed with dust lanes. Nonetheless, the 
photometry indicates that recent star formation has oc- 
curred in the SN environment. 

Figure 2 shows the color-color diagrams for the envi- 
ronment. The {U — B) vs. {B — V) diagram serves as 
the best diagnostic of the reddening and extinction in the 
SN's neighborhood, which, within the uncertainties in the 
stellar photometry, is Ay = 0-1.5 mag, with the mid- 
value Ay = 0.75 mag [E{B — V) ~ 0.24 mag, assuming a 
CardeUi, Clayton, & Mathis (1989) reddening law]. This 
is consistent with the locus of most of the stars on this 
diagram. We adopt this midvaluc, although the redden- 
ing in the environment appears to be variable and may 
at some locations exceed this value, up to Ay » 2 mag. 
Figure 3 shows the color-magnitude diagrams for the en- 
vironment. The brightest blue star on the diagrams is, of 
course, SN 1993J, which clearly has unusual colors, due 
to its emission-line spectrum. From these diagrams, Stars 
C and D are consistent with late 0-type to early B-type, 
likely on the main sequence (other possible fainter, young 
main-sequence stars arc also seen in the environment). 
Star A is quite red and is consistent with a K-type su- 
pergiant, while, from Figures 3a and 3b, Star B is either 
more evolved than Stars C and D or is experiencing greater 
reddening. The detected fainter stars, with V > 24.5 mag, 
tend to be red, owing to the relative lack of sensitivity of 
the U and B bands (see Table 1). 

The spectral energy distributions for Stars A-D are 
compared with that observed for the SN 1993 J progenitor 
in Figure 4. The archival images used by AHR were taken 
under a wide range of conditions, but, even for those ob- 
tained under the best conditions, the magnitude estimate 
for the progenitor still contained some or all of the light 
from these four stars. Star A contributed ~25% of the I- 
band flux measured by AHR for the progenitor, but a neg- 
ligible amount of blue light. The three other stars, B-D, 
contributed very little red light, but together contributed 
nearly 80% of the blue light and may be the source of the 
apparent excess flux at short wavelengths. One concern 



might be that background subtraction performed by AHR 
as part of the measurement process for the ground-based 
images already accounts for the mean contribution of these 
stars in an aperture, since the spatial distribution of stars 
within the aperture is similar to the distribution of stars 
in the overall background (due to the lack of a well-defined 
OB association around the progenitor). However, presum- 
ing that the mode is used to estimate the aperture photom- 
etry background in a sky annulus outside of the aperture, 
the fainter pixels in an image are given greater weight than 
the brighter ones concentrated in resolved stars. (Besides, 
AHR preferred PSF-fitting over aperture photometry in 
their measurements.) Therefore, the nearby stars would 
still contaminate the PSF of the progenitor star measured 
from the ground. 

To determine the likely energy distribution of the 
SN 1993J progenitor we must correct the AHR distribu- 
tion for this contamination, limits for which are shown in 
Figure 5. The sum of the light from the four stars is com- 
puted with (1) equal weight (dotted line) and (2) following 
a possibly more realistic case, in which the contaminating 
star contribution declines with angular distance from the 
progenitor (dot-dashed line). Because the central part of 
a ground-based PSF is similar to a Gaussian, the degree of 
contamination is likely to decline with distance in a similar 
manner. To weight the contribution from Stars A-D we 
use a Gaussian distribution, which is the sum of the unit- 
weight central distribution of the progenitor light and the 
distribution of the light from the contaminating star, 
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where Vi is the distance of the contaminating star from the 
progenitor and a is the width of the distribution [abbre- 
viated "E(F*g)" in Figure 5]. Since, for a Gaussian, the 
FWHM = 2.35(7, and the FWHM is a function of wave- 
length, a is also wavelength-dependent. From the AHR 
observational data, excluding the V and R measurements 
of Perelmuter (1993; AHR note that these measurements 
are several a above the mean in these bands), the values 
for (J, weighted by the standard deviations in the measure- 
ments, are au = 0'.'85, ae = 0'.'61, cry = 0'.'69, gr = 0'.'77, 
and (Tj = 0'.'53. The corrected progenitor energy distribu- 
tions are shown in Figure 5 as the long-dashed line (equal 
weight) and the short-dashed line (Gaussian weight). 

The actual energy distribution of the progenitor star 
probably falls somewhere within these limits. In the ul- 
traviolet (UV) the total flux from the surrounding stars 
is likely ~1(t from the AHR f-band flux. Thus, within 
the errors, the blue stars near the SN could account for all 
the UV flux of the progenitor measured by AHR, and the 
corrected [/-band flux from the progenitor is essentially 
undetermined. To possibly better estimate the progen- 
itor's energy distribution at redder bands, we simulated 
the light of the progenitor, plus the contamination from 
the nearby stars, as would be measured from the ground. 

We did this by first subtracting the SN PSF from the 
2001 HST F450W (B), F555W {V), and F814W (/) im- 
ages (when the SN was the faintest), employing a TinyTim 
model PSF (Krist 1995) in each band and the routines 

^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc. , under cooperative agreement with the National Science Foundation. 
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ALLSTAR and SUBSTAR (Stetson 1987, 1992) within 
IRAF/DAOPHOT. Wo then added at the SN position a 
fake star, representing the progenitor, using ADDSTAR 
within IRAF/DAOPHOT. We did this in a series of fake 
star trials, with a range of different input apparent mag- 
nitudes. These trial images were then convolved with a 
two-dimensional Gaussian with the wavelength-dependent 
values of a above, to simulate ground-based observing con- 
ditions. (These convolved images also then contained the 
light of the contaminating stars.) The model PSF was also 
convolved with the same Gaussian and subsequently ap- 
plied to each of the convolved trial images in ALLSTAR, 
to obtain magnitudes for each of the fake star trials in each 
band. We then determined which convolved fake progen- 
itor magnitude best matched that measured by AHR in 
each band for the actual progenitor. The input magni- 
tude which resulted in this best match likely provided the 
best estimate of the actual progenitor's brightness in each 
band. 

We show in Figures Ga 6c the input fake, or hypotheti- 
cal, progenitor images, before convolution, which resulted 
in the best matches to the AHR magnitudes after convolu- 
tion. [Note the presence of the light echo (Liu et al. 2002; 
Sugerman & Crotts 2002), particularly in the F450W and 
F555W images.] The resulting magnitude estimates for a 
hypothetical progenitor star are B « 22.4, V « 21.6, and 
1 K. 19.7 mag. In the F450W image the faint, very nearby 
Stars E and F are only ^1/4 the hypothetical progeni- 
tor star's brightness (some negative pixels resulting from 
oversubtraction of the SN are also seen) . The fake progen- 
itor in the / band has nearly the same brightness as the 
SN itself in the 2001 F814W image and is only ~1 mag 
brighter than Star A. In a sense, these final trial images 
simulate pre-SN WFPC2 images, had it been possible for 
them to have been obtained. 

We show in Figure 7a the final Gaussian-convolved fake 
progenitor image for F814W (/) as an example (compare 
this with Figure 6c). In general, the final test images in 
all three bands closely resemble the ground-based images 
shown by AHR (their Figure 1). Note the elliptical shape 
of the fake progenitor in the convolved /-band image, due 
primarily to the presence of Star A. (For the /-band trials. 
Star A dominated the resulting convolvcxi stellar profile 
for a faint fake progenitor, but the profile resolved into 
two stars as the trials were increased in brightness.) We 
show in Figure 7b the output image from ALLSTAR, with 
the fake progenitor subtracted. Star A becomes far more 
obvious after subtraction, although this star, as well as 
Stars B-D (which are much fainter at /), are all within 
the seeing disk of the convolved fake progenitor. (Compare 
Figure 7b with Figure 1 in AHR; their PSF-subtractcd /- 
band image, with the progenitor removed, clearly shows 
the presence of an object which we suspect is Star A.) 
Our simulation demonstrates that the ground-based mea- 
surements, even under good seeing conditions, contained 
not only the light of the progenitor star, but also that of 
the nearby contaminating stars. 

The formal uncertainties in the magnitude estimates for 
the progenitor, mostly due to aperture effects, small PSF 
mismatches, slight PSF oversubtraction, small disagree- 
ments between HSTphot and DAOPHOT results, and 
transformations to standard magnitudes, are estimated to 



be ^0.30 mag in each band. We have estimated that vari- 
able seeing in B and Y (corresponding to the range in 
seeing for the data in AHR) leads to an additional un- 
certainty of ~0.1 mag. Furthermore, the uncertainties in 
the V and, particularly, B magnitudes could be somewhat 
larger, since the faint stars E and F (and possibly others) 
very close to the progenitor are comparable to the progen- 
itor's brightness. Moreover, the relatively blue light echo 
(Liu et al. 2002; Sugerman & Crotts 2002), although faint 
and almost 2" from the SN, may contaminate the con- 
volved fake progenitor's brightness in B and V to some 
extent. The magnitude estimates are shown in Figures 3 
and 5. To within the uncertainties, the simulated magni- 
tudes agree fairly well with the limits shown in Figure 5, 
although the consistency is somewhat better in B and / 
than in V . 

Also shown in Figure 3 are solar-metallicity isochrones 
from Bertelh et al. (1994) for 8 Myr and 16 Myr, red- 
dened by E{B -V) = 0.24 mag and adjusted for the M81 
distance. The 16 Myr isochrone fits reasonably well the 
locus of Star A. If we assume that the SN progenitor and 
Star A were coeval, this age provides at least a lower limit 
to the progenitor's mass, as 16 Myr is roughly the main- 
sequence lifetime of an 12-14 Mq star (Schaller et al. 1992; 
Bertelli et al. 1994). With this assumption, we would ex- 
pect that the SN progenitor was at least as massive as 
Star A, since Star A is likely a red supergiant and there- 
fore also will likely explode as a supernova. The locus of 
the fake progenitor star on the diagrams is reasonably well 
fit, particidarly in the red, by the 8 Myr isochrone, which 
is roughly the main-sequence lifetime of a 21-24 Mq star. 
However, in the blue, the fake progenitor is somewhat bet- 
ter fit by the blue loop of the 16 Myr isochrones (again, 
the actual uncertainties in B and V for the fake progenitor 
are likely larger than the formal uncertainties). Bluer col- 
ors for the progenitor would agree with the evidence that 
the progenitor lost most of its hydrogen envelope (Filip- 
penko et al. 1994) and with a possible transition from red 
to blue supergiant during the late stages of the progenitor 
evolution (e.g., Immler, Aschenbach, & Wang 2001). 

Additionally, we can assign an approximate spectral 
type and luminosity class to the uncontaminated energy 
distribution for the progenitor, through construction of 
model atmospheres for supergiants using the ATLAS 12 
code (Kurucz 1993). One can see that the model colors 
for an early K-type (KO) supergiant agree reasonably well 
with the overall distribution. Because the progenitor lost 
most of its hydrogen envelope (e.g., Filippenko et al. 1994) 
and the progenitor atmosphere was helium-rich (e.g., Gar- 
navich & Ann 1994), we have calculated the KO I broad- 
band colors for three helium abimdanccs: solar, 75% he- 
lium, and 90% helium (by number). As expected, the col- 
ors for the solar-abundance model are very close to those 
given by Bessell (1990) from the Vilnius spectra. Unfortu- 
nately, the three models primarily differ in the UV, where 
the correctcid encirgy distribution of the progenitor is rela- 
tively undetermined. At longer wavelengths the colors of 
all three KO I star models are in good agreement with the 
probable energy distribution of the progenitor. 

Finally, after accounting for the contamination by 
nearby stars, we can place a limit on the brightness of any 
blue companion to the SN 1993J progenitor. The models 
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show that the KO I progenitor had U ~ 2A mag, thus con- 
tributing very httle to the flux measured by AHR in that 
band. A close hypothetical progenitor companion with 
U < 22.4 mag, added to the light from Stars A-D, would 
cause the total UV flux to exceed the AHR measurement 
for the progenitor by ~1(t. Conservatively, any progenitor 
companion could not have been brighter than U = 22.0 
mag and still be consistent with those measurements. The 
AHR uncertainty in the B band is small, so adding an ad- 
ditional star oi B < 23.0 mag to the ones accounted for 
here would exceed the observed flux by ~2(T for a KO I pro- 
genitor. Thus, the progenitor companion must have been 
less luminous than Mb ~ — 5.6 mag, assuming Ay = 0.75 
mag. Since 06 V and later-type dwarfs have Mb > —5.4 
mag (My > —5.0 mag; Humphreys & McElroy 1984) and 
masses of up to 30 Mq , we find that any hypothetical com- 
panion actually could have been of very high mass and still 
remained undetected. 

4. CONCLUSIONS 

We have used HST images of SN 1993J after outburst 
to search for stars within a ground-based resolution ele- 
ment, which may have contaminated the archival ground- 
based observations of the SN progenitor. We find four stars 
within 2'.'5 of the progenitor which were at least partially 
included in the progenitor energy distribution determined 
by AHR. Correcting for this contamination we find the 
following. 

1) The progenitor may have had apparent V — 21.6 ±0.3, 
B - V ^ 0.8 ± 0.6, and V" - / = 1.9 ± 0.5 mag, 
implying an absolute magnitude My ~ —7.0 ± 0.4 
for Av — 0.75 mag (the uncertainty in the abso- 
lute magnitude includes the estimated photomet- 
ric uncertainty, plus the uncertainty in the distance 
modulus). The brightness and colors are consis- 
tent with those of a reddened early K-type super- 
giant, although possibly somewhat too blue in the 
blue bands. The mass of the progenitor can be con- 
strained by the possible mass for the nearby red su- 
pergiant. Star A, i.e., '^13 Mq, and the age and mass 
of a hypothetical progenitor suggested by its bright- 
ness and colors, i.e., ~22 Mq. This range, 13-22 
Mq, is consistent with other estimates for the pro- 
genitor mass in the literature (e.g., Podsiadlowski 



et al. 1993; Woosley et al. 1994; Young, Baron, & 
Branch 1995; Iwamoto et al. 1997). 

2) Model atmospheres with a wide range of helium abun- 

dances fit the observed progenitor energy distribu- 
tion redward of 4000 A. Although at shorter wave- 
lengths the models show a greater spread in color 
with helium fraction, the data are insufficient to dif- 
ferentiate between them. 

3) The blue excess seen by AHR can be fully accounted 

for by the nearby stars. The progenitor was not a 
member of a significant OB association. Any com- 
panion to the progenitor must have had Mb > —5.6 
mag for reasonable assumptions of the reddening and 
distance to the SN. The mass of the companion is not 
very strongly constrained, but must be less than 30 
Mq . This is consistent with the binary mass-transfer 
scenario, even for conservative mass exchange (e.g., 
Podsiadlowski et al. 1993). 

SN 1993J appears to be interacting with circumstel- 
lar material ejected by the progenitor during its evolution 
(e.g., Filippenko et al. 1994; Van Dyk et al. 1994; Patat, 
Chugai, & Mazzah 1995; Matheson et al. 2000). The en- 
ergy from this interaction continues to power the optical 
light curve, which in 1995 slowed its decline to less than 
2 millimag d~^ (Garnavich et al. 1995). Comparing the 
magnitudes derived from the HST WFPC2 images in 1995 
and 2001, the SN is still declining at only '^0.2 mag yr"-'^. 
Recent ACS multi-band images have been obtained with 
HST by program GO-9353, but even with the superior 
sensitivity and resolution of ACS, the SN is likely still too 
bright to successfully isolate a companion or other stars in 
the immediate SN environment (<0'.'l-0'.'2). Only when 
the SN sufficiently fades, probably not for a few more HST 
cycles, can the high-resolution imaging be used to search 
for a companion to the SN progenitor and better determine 
the progenitor's nature. 
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Table 1 

Summary of HST Observations 



Date 


Band 


Exp. Time 




Program 


(UT) 




(s) 


(mag) 




1994 Apr 18 


F336W 


1200 


23.7 


GO-5480 




F439W 


600 


24.5 






F555W 


300 


25.5 






F675W 


300 


25.0 






F814W 


300 


24.5 




1995 Jan 31 


F255W 


2400 


22.1 


GO-6139 




F336W 


1160 


23.8 






F439W 


1200 


25.4 






F555W 


900 


26.6 






F675W 


900 


26.1 






F814W 


900 


25.5 




2001 Jun 4 


F450W 


2000 


26.7 


GO-9073 




F555W 


2000 


26.8 






F814W 


2000 


25.2 





'^These are the approximate HSTphot 3cr detection limits, in magnitudes, for the PC chip (programs GO-5480 
and GO-6139), and for the WF4 chip (program GO-9073). Magnitudes are given in the HST flight system. 



Table 2 

Measurements of Bright Stars Within 2.5" of SN 1993 J 



Star 




13 A a 
Jr^. A. 


Year 


u 


JD 


T/ 
V 




T 
1 




II 


o 




(F366W) 


(F439W)'' 


(F555W) 


(F675W) 


(F814W) 


A 


0.73 


348 


1994 


23.82(14) 


23.84(27) 


22.63(04) 


21.66(03) 


20.67(03) 








1995 




24.38(18) 


22.86(03) 


22.03(02) 


20.80(01) 








2001 




23.80(03) 


22.54(01) 




20.60(01) 








AVE. 


23.82(14): 


23.82(04) 


22.58(01) 


21.84(03) 


20.70(01) 


B 


1.44 


184 


1994 


23.01(14) 


23.12(11) 


22.92(05) 


23.05(07) 


23.16(14) 








1995 


22.94(14) 


23.27(08) 


23.02(03) 


23.12(04) 


23.06(05) 








2001 




23.17(02) 


23.12(02) 




22.84(10) 








AVE. 


22.98(14) 


23.17(03) 


23.07(03) 


23.10(05) 


23.03(07) 


C 


1.22 


353 


1994 


23.03(11) 


23.98(22) 


23.74(08) 


23.67(12) 


23.69(24) 








1995 


22.97(12) 


23.76(12) 


23.78(04) 


23.88(08) 


23.55(09) 








2001 




23.86(03) 


23.77(03) 




23.46(11) 








AVP 






ZO. ( / I I 


9*^ S9/'nQ'l 


9*^ i^*??! 1 \ 


D 


1.19 


149 


1994 


22.90(10) 


23.78(19) 


23.80(08) 


23.89(14) 


24.12(29) 








1995 


23.04(12) 


23.78(11) 


23.96(05) 


24.04(08) 


23.90(10) 








2001 




23.78(03) 


23.90(03) 




24.19(14) 








AVE. 


22.96(11) 


23.78(04) 


23.90(04) 


24.00(09) 


24.01(13) 


E 


0.32 


335 


1994 


















1995 




23.84(14) 


24.04(10) 


23.68(12) 


23.82(14) 








2001 




24.59(23) 


24.32(07) 






F 


0.29 


282 


1994 


















1995 


24.35(34) 


24.03(13) 


24.87(09) 


23.77(07) 


23.96(11) 








2001 













Note: The values in parentheses are the la uncertainties in the last two digits of the magnitudes. 
^Position angle, degrees from north through east. 

''The 2001 observations by GO-9073 were made through the F450W, rather than the F439W, filter. 



Table 3 

HST Photometry of SN 199:^.1 



Date 


UV 


U 


B 


V 


R 


/ 


(UT) 


(F255W) 


(F366W) 


(F439W)^ 


(F555W) 


(F675W) 


(F814W) 


1994 Apr 18 




17.98 


18.22 


17.78 


16.99 


17.21 


1995 Jan 31 


17.30 


18.55 


19.10 


18.75 


17.90 


18.63 


2001 .Tuu 04 






20.4<S 


20. .32 




19.75 



Note: The la imcertainties in the magnitudes are all <0.01 mag. 

^The 2001 observations by GO-9073 were made through the F450W, rather than the F439W, filter. 
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Fig. 1.— The sum of the HST images obtained in 1995 through the F336W, F439W, F555W, F675W, and F814W filters, 

showing the environment around SN 1993J in M81. The circle centered on the supernova has a radius of 2'.'5, which 
is approximately the full-width near zero intensity (FWZI) of a Gaussian seeing disk with full-width at half maximum 
(FWHM) = 1'.'5 (typical seeing for observations analyzed by AHR). The six brightest stars in the SN environment are 
labelled. The four brightest, A-D, with V < 24 mag, most likely contaminated the ground-based measurements of the 
SN progenitor by AHR. 




-0.5 0.5 1 1.5 

B-V (mag) 

Fig. 2. — The color-color diagrams for the SN 1993J environment (an area 2'.'5 in radius, centered on the SN), measured 
from the 1995 HST images, (a) The {U — B, B — V) diagram; and (b) the {V — R, V — I) diagram. In both diagrams. 
Stars A-D are represented by open triangles and labelled, whereas open squares represent other stars in the environment. 
The SN {five-pointed star) is also labelled. In (a), Star A is not shown, due to the lack of a [/-band measurement. In 
both diagrams the main-sequence locus is shown {solid line), as well as the reddening vector {dashed line), assuming a 
Cardelli et al. (1989) reddening law. 
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Fig. 2. — (Continued.) 
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18 



20 



<20 22 

CO 

s 



24 - 



26 - 



T 1 1 1 1 1 1 1 r 



1 1 1 \ \ 1 1 1 1 1 r 

SN 1993J environment 

(a) 





16 Myr, 
\ 



II ^rm^ 




Av=0 mag 




Av=l mag 



1 



1 



1 
U-B (mag) 



Fig. 3. — The color-magnitude diagrams for the SN 1993 J environment (an area ~ 2'.'5 radius, centered on the SN), 
measured from the 1995 HST images, (a) The {U - B, V) diagram; (h) the {B - V, V) diagram; and (c) the {V - I, V) 
diagram. In all three diagrams, the five-pointed star with V = 18.75 mag is the SN, and Stars A D are represented by 
open triangles and labelled, whereas open squares represent other stars in the environment. (In (a) and (b), these four 
other blue stars are to the north and northwest of the SN. In (c), the several faint, red stars are distributed throughout 
the environment. The faintest detected stars are the reddest, due to the higher sensitivity of the redder 1995 images.) 
The main sequence locus is the solid track; an 8 Myr isochrone with solar metallicity from Bertelli et al. (1994), reddened 
assuming Ay = 0.75 mag and adjusted for the M81 distance from Freedman et al. (1994), is shown as the long-dashed 
track, and a 16 Myr isochrone is shown as the short-dashed track. A hypothetical progenitor star is represented by the 
open circle. The reddening vector, following a Cardelli et al. (1989) reddening law, is also shown. 



11 




Fig. 3. — (Continued.) 



18 
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20 



SN 1993J environment 

(c) 



> 



cio 22 

CO 



24 - 



I 



/ 



I 



^ 8 Myr 



16 Myr 



Ay=0 mag 



. 1 min ^I I I 



'k^=\ mag 



26 - 



£3= 



t3- 



H - 



J L 



J I L 



J I L 



J I \ I I I I I I I I I L 







0.5 1 1.5 

V-I (mag) 



2.5 



3. — (Continued.) 
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Fig. 4. — The energy distributions of the four bright stars, A D, within 2'.'5 of SN 1993J, compared to the energy 
distribution for the SN progenitor from AHR. The short-dashed line shows the sum of all four stars. 
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20 



(0 



24 



"1 — 1 — I — I — I — I — I — I — I — I — i — I — 1 — I — I — I — 1 — I — I — I — r 



sum 



.^2(F*g) 




AHR-sum 

AHR-i:(F*g) 

AKO I solar, Av=0.75 mag 
■ KO I 0.75 He, Av=0.75 mag 
♦KO I 0.90 He. Av=0.75 mag 



I I I I I I I I I I I I I I I I I I I 



J I I I I L 



4000 5000 6000 7000 8000 

Wavelength (A) 

Fig. 5. — The progenitor energy distribution corrected for the contamination by the four nearby stars. The solid line 
shows the progenitor briglitness from AHR, using archival ground-based observations; the dotted line represents the equal- 
weighted sum of the light from all four stars; the dot-dashed line represents E(F*g), the sum of the flux, F, from each 
of the four stars, weighted by their distance from the SN, assuming a one-dimensional Gaussian distribution, g, with 
wavelength-dependent width, a (see text). The long-dashed line represents the progenitor energy distribution corrected 
by the equal-weight sum, while the short-dashed line represents the distribution corrected by the Gaussian-weighted sum. 
The former is likely a lower limit to the true progenitor spectrum, while the latter is likely an upper limit. The open 
squares represent the light from a fake, hypothetical progenitor (see text). The solid symbols show the energy distribution 
for model K-type supergiant atmospheres with a range of possible helium abundances, reddened assuming Ay = 0.75 
mag. 
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Fig. 6.— (a) The F450W (B) HST image from 2001 June of the SN 1993J environment, after subtraction of the SN 
via fitting of a TinyTim PSF (Krist 1995) within DAOPHOT/ALLSTAR and the addition of a fake progenitor star with 
B Ki 22 A mag. (h) The F555W {V) HST image from 2001 June, after SN subtraction and the addition of a fake progenitor 
star with V « 21.6 mag. (c) The F814W (/) HST image from 2001 June, after SN subtraction and the addition of a fake 
progenitor star with I w 19.7 mag. The light echo (Liu et al. 2002; Sugerman & Crotts 2002) is indicated in (a) and (b) 
by an arrow. 



Fig. 7. — (a) The F814W (/) HST image from 2001 June with the fake progenitor star (see Figure 6), after convolution 
with a two-dimensional Gaussian with a = 0'.'53. (b) The same as (a), except after subtraction of the fake convolved 
progenitor via fitting of the TinyTim PSF, also convolved with the same Gaussian, within ALLSTAR. The position of 
Star A is indicated. 
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